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^ ; Abstract 

■ Three ways of constructing a non-Hermitian matrix with possible all real eigenvalues are dis- 
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cussed. They are VT symmetry, pseudo-Hermiticity, and generalized VT symmetry. Parameter 
counting is provided for each class. All three classes of matrices have more real parameters than 
a Hermitian matrix with the same dimension. The generalized VT -symmetric matrices are most 
general among the three. All self-adjoint matrices process a generalized VT symmetry. For a given 
^ matrix, it can be both "PT-symmetric and 'P-pseudo-Hermitian with respect to some V operators. 

cr 

The relation between corresponding V and V operators is established. The Jordan block structures 
of each class are discussed. Explicit examples in 2 x 2 are shown. 

T— I ! 

oo 

7— I 

<N 
(N 



I. INTRODUCTION 

Since Bender and Boettcher published their seminal paper in 1998 jlj], the field of so-called 
PT-symmetric quantum mechanics grows rapidly. Two review articles have been published 

as. 

Applications in optics are widely studied in both theoretical predictions flfl and 
experimental realizations |s-12|. 



The essential idea of this new formulation of quantum mechanics is to allow the Hamil- 
tonian to be apparently non-Hermitian, yet the entire spectrum is still real. In finite di- 
mensions, this means that we may represent the Hamiltonian of a non-dissipative quantum 
system by a non-Hermitian matrix. Namely, a matrix H with all real eigenvalues but 
H ^ . In this paper, we try to classify all non-Hermitian matrices with possible all real 
eigenvalues. We hope our results will provide a reference on 'PT-symmetric matrices to the 
community. This work can be considered as a extension with more general parity operators 
studied in Ref. 13] and in higher dimensions discussed in Ref. \\\ . 

For a diagonalizable matrix H to have all real eigenvalues, a necessary and sufficient 
condition is that there exist a Hermitian matrix W = such that 

WE = H^W (1) 

n 

and all the eigenvalues of W are positive definite |l5|. If H represents the Hamiltonian of 
a quantum system, then W is the metric operator which defines the inner-product of the 
corresponding Hilbert space jl^ . 

(•,-) = <-im (2) 

where (-| = Eq. ([1]) immediately leads to 

(;H-) = (H;-). (3) 

That is, H is self-adjoint. In the conventional quantum mechanics, the Hamiltonian can 
always be represented by a Hermitian matrix. In this case, on can always choose W to be 
the unity matrix. And the inner-product reduces to the familiar Dirac's one between bras 
and kets. 

In the first paper on "PT-symmetric quantum mechanics, Bender and Boettcher identified 
the combined symmetry of P(parity) and 7~(time reversal) as the generalization of Hermitic- 
ity jl|. They coined the term "PT-symmetric Hamiltonian" for a Hamiltonian commuting 
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with PT, 

[VT,H]=0. (4) 

Later, Bender et al realized that the metric operator in PT-symmetric quantum mechanics 
is non-trivial 17|. To show the unitary time evolution, they introduced a so-called "C 
operator' 1 and defined a CPT inner-product. Translating to the current terminology, Bender 
et al defined the metric operator as W = VC. 

However, the PT symmetry does not guarantee the Hamiltonian to have an entirely real 
spectrum. It only ensures a real secular equation of H because 

(PT) 2 = 1. (5) 

This means that all the eigenvalues of H either are real or form complex conjugate pairs. 
If all the eigenvalues are real, we call the PT symmetry is unbroken. Otherwise, the PT 
symmetry is spontaneously broken. Not all PT-symmetric Hamiltonians are diagonalizable 
either. The Jordan block structure may be found if there is degeneracy. 

In Sec. m we study the PT-symmetric matrices by carefully defining the parity operator. 



Comparing to an early study in Ref . [13|] , our formulation here is more general and it allows 
us to study the non-symmetric parity and Hamiltonian operators. We will show that an 
(m+n) x (m+n) PT-symmetric matrix has (m+n) 2 +2mn real parameters. For comparison, 
an (m + n) x (m + n) Hermitian matrix has (m + n) 2 real parameters. We will also show 
that all possible Jordan block structure may be found in the PT-symmetric matrices. 

If H has complex eigenvalues or H is not diagonalizable, the metric W in ([1]) stops 
being well defined: Either W becomes non-Hermitian, or the eigenvalues of W are no longer 
positive-definite. If one maintains the Hermiticity of W but allows it to have negative 
eigenvalues, the linear space equipped with the indefinite metric is called the Pontrjagin 
space (a finite-dimensional version of the Krein space). In the literature, the Hamiltonian 
H with an indefinite-metric is called pseudo-Hermitian j^]. 

In Sec. IHIt we discuss the P-pseudo-Hermitian matrices where P plays the role of an 
indefinite metric. This section can be considered as an extension of the early study in 
Ref. [l^l to higher dimensions. We show that an (m + n) x (m + ri) P-pseudo-Hermitian 
matrix has (m + n) 2 + Iran real parameters, same as in a PT-symmetric matrix in the 
same dimension. A Hermitian matrix can always be considered as a special case of P- 
pseudo-Hermitian matrices. And all possible Jordan block structures may be found in the 



P-pseudo-Hermitian matrices. 

It is commonly believed that the VT symmetry is neither necessary nor sufficient condi- 
tion for a Hamiltonian to have an entire real spectrum. This is not true if we generalize the 
VT symmetry to be a combined involutory operation VT. In Sec. \W\ we show that these 
generalized PT-symmetric matrices contain all self-adjoint matrices. The generalized VT 
symmetry is the necessary condition for a matrix to have all real eigenvalues. In the case 
of N x N, a generalized PT-symmetric matrix has 2N 2 — N real parameters, same as in 
an N x N self-adjoint matrix. All the PT-symmetric matrices studied in Sec. [IT] are special 
cases of generalized PT-symmetric matrices. If one wants to model a finite-dimensional non- 
dissipative quantum system in the most general way, he/she should start with a generalized 
PT-symmetric matrix Hamiltonian. 

The explicit 2x2 matrices are presented for each class. In this dimension, it is interesting 
to know that all three classes coincide: A PT-symmetric matrix is also P-pseudo-Hermitian, 
and it exhibits the generalize VT symmetry as well. Only in higher dimensions, three classes 
may contain different matrices. And both PT-symmetric and P-pseudo-Hermitian matrices 
are special cases of the generalized PT-symmetric ones. 

In the concluding Sec. |Vj we summarize the parameter counting for different types of 
matrices in Table [TJ 

II. VT SYMMETRY 

In this section, we study PT-symmetric matrices. Let us start with the definitions of the 
parity operator P and the time reversal operator T. Since the time reversal operator is an 
involution and anti-linear, we define it simply as the complex conjugation, 



That is, for a matrix (operator), TAT = A*, for a vector (state), T\ip) = |"0)*> where * 
represents complex conjugation. 

The parity is a linear operator. In finite dimensions, we define it as a matrix, P. As 
usual, we demands that the time reversal and the parity commute from each other, 



T = * T 2 = 1. 



(6) 



[p,n = o. 



(7) 
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This commutation relation immediately leads to that V must be real. Together with the 
involutory property, V satisfies 

V = V* and V 2 = 1. (8) 



Our definition of the parity operator is more general than the one studied in Ref. [13|, in 
which the parity operator is always symmetric because only (complex) symmetric matrix 
Hamiltonians were considered. Here, we keep V to be real but allow it take non-symmetric 
form. This generalization gives us more parameters in both the parity operator and the 
Hamiltonian. 

We say that a matrix Hamiltonian H is T-'T-symmetric if it commutes with the combi- 
nation of VT. With our definitions of V and T here, the PT-symmetric condition is 

[PT, H} = VH = H*V. (9) 

Because both V and T are involutions and they commute, the combination of VT is also an 
involution, {VT) 2 = 1. One consequence of this result is that the eigenvalues of H can either 
be real or form complex conjugate pairs. The VT symmetry is not a sufficient condition for 
a matrix to process only real eigenvalues. If a PT-symmetric matrix Hamiltonian H has 
an entirely real spectrum, we call that the VT symmetry is unbroken. In this case, all the 
eigenstates of H are also eigenstates of VT. 

H\E n ) = E n \E n ), VT\E n )=V\E n y = \ n \E n ) (10) 

with E n to be real E n = E* and A n to have modulus one, A n = e ia ". One can always choose 
the phase of the eigenvectors such that A n = 1. When complex eigenvalues appear in H, 
even though H is still "PT-symmetric, the eigenstates of H are no longer eigenstates of VT. 
In this sense, we call that the VT symmetry is spontaneously broken. 



A. Similarity transformations 

Because V is an involution, its eigenvalues can take only the values of ±1. We may 
use the number of positive and negative eigenvalues to classify a finite-dimensional parity 
operator. Let us denote an m + n-dimensional parity operator with m positive eigenvalues 
and n negative eigenvalues as V(m,n). All V(m,n) operators with the same m and n can 
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be linked by a similarity transformation. Since P is a real matrix, we consider only the real 
transformations. All these transformations form a general linear group, 

V(m,n) = RVoim^n)^ 1 , with R e GL(m + n, E). (11) 

This transformation maintains both reality and involutory properties of "P. There are (m + 
n) 2 real parameters in a transformation matrix i?, but not all of them enter V as we will 
see later. 

The same transformation preserves VT symmetry. That is, if Hq is Po7~-symmetric, then 

H = RHqR- 1 (12) 

is PT-symmetric: 

V H = H*V VH = H*V. (13) 

Just as shown in Ref. 13j, we may use this property to construct PT-symmetric matrices. 



Without loss of generality, we choose the parity operator P to be diagonal, 

V (m,n) = Diag{l, ---,1, -1, ■■•,-!} = . (14) 

y -InxnJ 

It is easy to find that for an (m + n) x (m + n) matrix H Q to be PoT- symmetric, it must 
have the block form 

„ / ^mxm iB mxn \ 

H = , (15) 

V iC n xni Dnxn I 

where A, B, C, and D are all real matrices. There are m 2 + n 2 + 2mn = (m + ra) 2 real 
parameters in H . 

To count the number of real parameters in a generic (m, n) parity operator, let us take a 
closer look at the similarity transformation in ( fTTj) . Many members of GL(m + n,M) leave 
Vo{m,n) invariant. By group theory, all of such transformations form a group, it is called 
the stabilizer subgroup or the little group. Suppose that R' is a member of this subgroup, 

R'V (m, n)^- 1 = P (m, n), & [#, P (m, n)] = 0. (16) 

It is easy to see that R' must be block diagonal, 

I Xmxm ° 
Y nx 



with X and Y to be real invertible matrices. Hence, the little group of Vo(m, n) is the direct 
sum of GL(ra, R) and GL(n, R). All the nontrivial similarity transformations which defines 
V belong to the coset of 

(18) 



GL(m + n, . 



GL(m,R) x GL(n, R) ' 

Thus, among (m + n) 2 real parameters in GL(m + n, R) group, m 2 in GL(m,R) subgroup 
and n 2 in GL(n, R) subgroup do not enter V. There are (m + n) 2 — m 2 — n 2 = 2mn real 
parameters in V. Combining with the parameters in H Q , we conclude that a PT-symmetric 
Hamiltonian has (m + n) 2 + 2mn real parameters. 

The metric operator can be found by solving the self-adjoint condition WqH = HqWq. 
When the parity is transformed to V — RVoR" 1 and the Hamiltonian changes to H = 
RHqR -1 , the metric operator transforms accordingly, 



W 



R-^WqR- 1 . 



(19) 



Obviously, the transformation preserves the Hermiticity of the metric operator. That is, 
W = if Wq = Wq. Although the eigenvalues of W may change under such a transfor- 
mation, but they remain to be positive. 

B. Jordan blocks 

Not all "PT-symmetric matrices are diagonalizable. They may form Jordan blocks. Let 
us denote J n (X) as the n-dimensional Jordan block with eigenvalue A, 

f \ 1 ••■ 
A ••■ 



J n (A) = Al nxn + N n 



••• A 1 
•■• A 



(20) 



/ nxn 



with N n the nilpotent matrix with units on the first upper diagonal and (A^)™ -1 ^ 0, 
(N n ) n = 0. For an A-dimensional matrix H, it is always similar to a direct sum of a series 
of Jordan blocks by a similarity transformation A, 



AHA' 1 = J ni (Ai) © J„ 2 (A 2 



with N = ri\ + ri2 + 



(21) 
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Only for the special case of n\ — rii = • • • = 1, H is diagonalizable. 

Since similarity transformations do not effect the diagonalizability, we only need to study 
the possible Jordan block structure in H . By choosing parameters in H Q , all possible Jordan 
blocks emerge in PT-symmetric matrices. For example, if we choose 

A = J m (A), 

••■ 

1 ••• 



B 



\ 

C = 0, 
D = J n (A), 



/ raxn 



then Hq is similar to J m+n (X): 



with 



AH A 1 = J m+n (X), 



A 



(22) 



(23) 



(24) 



Hnxn 

Similarly, all other possible Jordan block structures in (121j) can be constructed by properly 
choosing the real matrices A, B, C, and D in Ho. Of course, the construction here is only 
a sufficient condition to form Jordan blocks, it is not necessary. 



C. 2x2 



Now let us show the explicit matrices in the case of 2 x 2. We start with a simple diagonal 
parity operator, 

1 
-1 

where 03 is the third Pauli matrix. A PoT-symmetric matrix Hamiltonian can be parame- 
terized as 



(25) 



H = e<Jo + (ip, 7 sin 5, 7 cos 5) ■ cr 

e + 7cos5 — i(7sin5 — p) 
i(7sin5 + p) e — 7 cos 5 



(26) 



with four real parameters: e, 7, p, and S, where cr = (o"i, 02, 03) are Pauli matrices and 
0o = f-2x2- The eigenvalues of H are 



E± = e ± V7 2 - p 2 . (27) 

They are real when 7 2 > p 2 . It is the condition for VT symmetry not broken. The 
corresponding eigenvectors are 

^J^-^^ZE], (28) 
y 7e ld + ip =F V 7 2 - P 2 J 

where n± is the normalization constant. (The little bit complicated choice of eigenvectors 
here is to avoid the accidental degeneracy when p = ±7 sin 5.) It can be checked that \E±) 
are also eigenvectors of VqT with the eigenvalues to have modulus one. 

V T\E ± )=e ia ±\E±). (29) 

We can always choose the phases of the eigenvectors such that the eigenvalues of VqT to be 
one. 

Solving the self-adjoint condition in Eq. (JTJ), we get the form of the metric operator, 

Wo = u [700 + (0, v sin 5 — p cos 5, v cos 5 + p sin 5) ■ a] 

(7 + (v cos 6 + p sin 6) —i(v sin 5 — p cos 5) \ 
(30) 
i(v sin 5 — p cos 5) 7 — [v cos 5 + p sin 5) I 

where u and v are arbitrary real constants with the constraints w~y > and v 2 < 7 2 — p 2 . 
The eigenvalues of Wo are 

u± = u (7 ± Vp 2 + v 2 ) . (31) 
They are positive definite. With this metric operator, the eigenvectors of H are orthogonal 

(E+\W \E-) = = ( J EL|W / |£+) (32) 

and normalized to 

A4 = (E±\W \E±) = 4|n ± | 2 n 7v / 7 2 1 P 2 (V7 2 - P 2 ± «) • (33) 

Given it and one can always normalize the eigenvector by choosing proper n±. It is also 
true in the reversed way, for arbitrary non- vanishing n±, one can always tune u and v in the 
metric to normalized the eigenvectors. 
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GL(2,K) , . , , 

A generic matrix m the coset of — — — — — — — — can be parameterized either by 

GrL(l,M) x GL(l,ffi) 

R 1 (6 1 <p) = e-^ 3/2 e- Wa2/2 J (34) 

or by 

R 2 (9,<p) = e -^/2 e -e CT1 /2_ (35) 

Note that they are not equivalent. According to these transformations, we have the general 
V has the form either as 

/ cos6> e _¥, sm6» \ 
Vi = Ri(e,<p)V Rt\e,<p)= \ , (36) 

V e^ sin 9 — cos 9 I 

or 

/ cosh 9 e~ v sinh 9 \ 
V 2 = R 2 (9,<p)V R 2 1 (9^)= \ . (37) 

y -e v sinh 9 - cosh 9 J 
Two types of T-'T-symmetric Hamiltonians can be obtained by the corresponding transfor- 
mations, 

H, = R 1 (9,i f )H R^ 1 (9,i f ) 

e + 7 cos 5 cos 9 — ip sin 9 (7 cos 5 sin 9 — iy sin 5 + ip cos 9)e~ ip 
(7 cos 5 sin 9 + iy sin 5 + ip cos 9)e v e — 7 cos S cos 9 + ip sin 9 

- R 2 (9 1 <p)H R 2 - 1 (9 1 <p) 

e + 7 cos(5 + i9) -i[7sin(5 + i9) - pje"^ 
i[7sin(5 + id) + p]e v e - 7 cos(<5 + i9) 

Finally, the metric operators for each H transformed accordingly 

Wi = (Ri i y(9, ( p)W RT 1 (9,<f) 

[7 + cos 9(p sin 5 + v cos 5)} e v sin 9(p sin 5 + v cos 5) + i(p cos 5 — v sin 5) 

sin #(p sin 5 + v cos 5) — i(p cos 5 — v sin 5) [7 — cos 9(p sin 5 + i> cos 5)] e~ v 

(40) 

W 2 = {R 2 1 ) j (9,^W R 2 1 (9^) 

([7 cosh 6> + (p sin 5 + v cos 5)] e^ 7 sinh 9 + i(p cos 5 — v sin 5) \ 
(41) 
7 sinh — i(p cos 5 — v sin 5) [7 cosh 9 — (p sin 5 + w cos 5)] e ^ / 

Now let us see what goes wrong when a Jordan block is formed. Remember that similarity 
transformations do not change the diagonalizability of a matrix. Let us work in Vq frame. 



,(38) 



(39) 



a 



The matrix Ho in Eq. (126]) forms a Jordan block if and only if 7 2 = p 2 ^ 0. Because 
v 2 < 7 2 — p 2 , to avoid double-limit, we fix v — 0. Let us define a small positive parameter e 
by 

7 2 = p 2 (l- e ). (42) 

For arbitrary value of u, when e is small, the larger eigenvalue of Wo, oj > behaves well but 
the smaller one, u < is vanishing, 

u > ~ 2^7, and u;< ~ |u7e. (43) 

A similar ill behavior can be found in the normalization of eigenvectors of Ho, 

A4 ~ 4|n ± | 2 w 3 e. (44) 

When Ho has a Jordan block structure, one can only find one eigenvector. Another vector 
to span the linear space lives in the Jordan chain. The non-diagonalizable H has the form 

(e + 7cos<5 — sin 5) \ 

>K ' , with 7^0. (45) 

i7(l + sin5) e — 7 cos 5 I 

The eigenvalue of H Q is e and the only eigenvector is 

(1 — sin 5 \ 
(46) 
i cos 5 I 

where no is the normalization constant. It is obvious that |$o) is an eigenvector of VT, 

m* 

VT\%) = ^|$ ). (47) 
n 

The second vector in the Jordan chain satisfies 

(H -ea )\$ 1 ) = \$ ). (48) 

The solution is 

( l-sin<5 \ 
^ +a\%) (49) 

with a to be an arbitrary constant. When a is real, |$i) is also an eigenvector of VT with 
the same eigenvalue as | < &o}> 

n* 

VT\$i) = (50) 
n 
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III. PSEUDO-HERMITICITY 



In this section, we explore the so called pseudo-Hermitian matrices, an alternative for- 
mulation of non-Hermitian matrices with possibly all real eigenvalues. Let us start with a 
Hermitian involutory operator V: 

V = V\ and V 2 = 1. (51) 

A P-pseudo-Hermitian matrix Hamiltonian H satisfies 

VHV = H ] & VH = H ] P. (52) 

Clearly, in this formulation, the V operator plays the role of the metric operator except 
that V may have negative eigenvalues. In general, the linear vector space with an inner- 
product defined by an indefinite metric operator is called Krein space. In the case of finite 
dimensions, it is a Pontrjagin space. 



A. Unitary transformations 

A unitary transformation preserves the property for the V operator defined in Eq. (15T1) . 
Just like the method used in the previous section, we may start with a simple operator Vq. 
We then use unitary transformations to obtain all possible V operator. For the case of m + n 
dimensions, 

V(m, n) = UV (m,n)U\ with UeU(m + n). (53) 

There are (m + n) 2 real parameters in U. Again, we will see later that not all parameters 
enter V(m, n). 

If we can find a Po-pseudo-Hermitian matrix Hamiltonian Hq, then all the P-pseudo- 
Hermitian matrix Hamiltonians H can be obtained by the same unitary transformation, 

H = UH Ul (54) 

Without loss of generality, let us choose the Vo(m,n) operator to be the same as the 
parity operator in the previous section, 



Po(m,n) = Diag{l,-- -,1,-1, ■■•,-!}= mXm . (55) 
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A Po-pseudo-Hermitian matrix Hamiltonian has the block form 



#0=1 | , (56) 

i(-SmXn) D nxn 

with B to be arbitrary and A and D to be Hermitian: A = A> and D = D\ Obviously, Hq 
has m 2 + n 2 + 2mn = (m + n) 2 real parameters. 

Now, we perform the unitary transformation in (1531) to get the general V operator. Note 
that the little group of Vo{m, n) in U(m + n) will leave Vo{m, n) invariant. This little group 
is the direct sum of U(m) and U(n). Thus, among (m + n) 2 real parameters in the group of 
\J(m+n), m 2 in the subgroup of U(m) and n 2 in the subgroup of U(n) do not enter V. There 
are (m+n) 2 — m 2 — n 2 = 2mn real parameters in V . The coset space which parameterizes the 
nontrivial similarity transformations of the V operators is a so called complex Grassmann 
manifold, or complex Grassmannian, 

TT y ( ? + T ?,\ = G m (C m+ ") (57) 
U(m) x U(n) 

As coset space it can be parameterized by the Lie algebra elements of the type 









a = 




, 4eC mx "| 




\-» 1 °J 





u(m + n) Q (u(m) © u(n)) = < a G u(m + n) a = \ , b G C mxn > , (58) 

i.e. by ant i- Hermitian matrices a = — G C (m+n)x(m+n) with the u(m) and u(n) blocks 
removed from the block-diagonal. For the corresponding coset space elements U of the 
group it holds the representation 

, cos(vWx) b s[n ^ bx) 
U = e ax = m , xeR. 59 



The general P-pseudo-Hermitian Hamiltonian can be defined as 

H = UH UK (60) 

It has (m + n) 2 + 2mn real parameters. 

In general Hq is not Hermitian because the factor of i in front of B and B\ To make 
it Hermitian, one has to choose B — 0. The Hermitian limit of Hq has the block-diagonal 
form 

i A m xm 

&o = > ( 61 ) 

D r 
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where A = A\ D = . It has m 2 + n 2 real parameters. Using the unitary transformation 
in Eq. fl57|) . we get the general Hermitian matrix, 

h = Uh Ul (62) 

Now, the question is whether all (m + n) x (m + n) Hermitian matrices belong to h. They 
do. Here is the reason. Since A and D are Hermitian, they can be written as unitary 
transformations of diagonal matrices, 

A = U m a m U ] m with U m G U(m), 
D = U n d n Ul with U n e U(n), 

(63) 

where a m and c? n are diagonal matrices. The combination of U, U m , and U n re-ensemble 
all the members in U(m + n). Hence, h includes all Hermitian matrices with the same 
dimension. 

One advantage for P-pseudo-Hermiticity is that the V operator defines a Krein space 
spanned by the eigenvectors of H. 

(^,^ = (^|P|0). (64) 

All eigenvectors of H with different eigenvalues are orthogonal with respect to this inner 
product. However, the norm is not positive definite. Only if the P-pseudo-Hermitian matrix 
Hamiltonian H has an entire real spectrum and it is diagonalizable, a proper Hilbert space 
can be defined by a positive definite metric operator W as in Eq. (E]). Not only different 
eigenvectors of H are orthogonal respect to this W- inner-product, but also the eigenvectors 
can be normalized properly. In general, W is dynamic, namely if-dependent: WH = H^W. 
W transforms the same way when one applies the unitary transformation on a "P-pseudo- 
Hermitian matrix Hamiltonian H: 

W = UW Ul (65) 

B. Jordan blocks 

Not all "P-pseudo-Hermitian matrices with real eigenvalues are diagonalizable. They may 
form Jordan blocks. In iV x N space, one can always find a unitary transformation U such 
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that 



v = uv u j = & © Si 



with i + j + 



N, 



(66) 



where S n is the so called standard invohitory permutation (SIP), 



/ 



■■■ 



••• 1 

1 ■■• 

1 ••■ 



^ 2 = 11 



(67) 



with units on the skew-diagonal. They provide the intertwining metrics for Jordan blocks 
in Eq. (P20D. For A € R, 



SrJnWSn = 4(A)- 



(68) 



Obviously, it holds that 

Tr(S 2n ) = 0, Tr(S 2n+1 ) = 1. (69) 

So that (because of S% = l nX n) we should have the similarity relations S^n — Vo(n,n), 
S2n+i — Vo(n + l,n). [This follows from the invariance of the trace under similarity trans- 
formations and the fact that Tr(Po(^; n )) = 0, Tr("Po( n + l? 7 ^)) = 1-] A simple calculation 
gives 



S2n 

Q 



qV (n,n)q \ 

1 / lnxn S n 
\ S n l nxn 



M/4 b 



o -s n 
s n 



1 / lnxn <S( 

I -S n 1 



(70) 



n -"-nxn 



as well as 



SWh-i = QVo(n + l,n)g 1 , 



1 

7! 



/ n n — <? ^ 



V2 

5n lnxn 



>/4 ft 



^ -S n ^ 



s n 
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V2 



-U-nxn u J n 



\/2 

-S n Inxn 



V 



(71) 



Here we have used the fact that 



S n 

Sn 

-1 

1 



-1 

1 



-1 

1 



L 2x2; 



(72) 



, - 1 

i.e. that | | can be considered as 2-dimensional realification of the imaginary unit i 

1 

(a representation of i in terms of a real 2x2 matrix). In other words, we can identify 



and calculate 



-1 

1 



cos(x)]l nxri + isin(a;)5' n 



what for x = 7r/4 and via Eq. (173]) just gives Eq. flTTj) . 

It can be verified that the following matrix Hamiltonian is P-pseudo-Hermitian, 



H = UH U j = Hi® Hi 



with i + j 



where H n has a structure looks like a "rotated Hermitian matrix". 



H„ 



an + ifen oi2 + ifoi2 

a 2 l + i&21 «22 + 



' «l,n-l + l&l,n-l a ln 

^2,71-1 Ol,n-l — i&l,n-l 



0-22 — 1^22 °<12 — 1^12 

a 2 i - ifoi a n - ib n J 



(73) 



(74) 



(75) 



(76) 



3-71-1,1 + i^n-1,1 a n-l,2 
y a n l On-1,1 — i^n-1,1 ' 

with all djj and bij to be real parameters. 

The advantage for this parametrization is that the Jordan block can easily be constructed. 
If we set an = a 2 2 = ■ ■ ■ = A, a\ 2 = a 2 % = ■•• = !, and all others to be zero, H n becomes 
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a Jordan block of order n. Thus, we have shown that all possible Jordan blocks emerge in 
P-pseudo-Hermitian matrices. Of course, the construction here is only a sufficient condition 
to form Jordan blocks, it is not necessary. 

When Jordan block is formed, H is no longer diagonalizable and the metric operator 
cannot be properly defined. We will show explicit examples in 2 x 2 later. 



C. Relation between VT symmetry and pseudo-Hermiticity 

In the earlier papers about PT-symmetric matrix Hamiltonians, only real-symmetric 
parity operators are considered 0, Q, H] • In this case, the parity operator in VT symme- 
try is the same as the V operator in P-pseudo-Hermiticity. The transformations between 
different real-symmetric V operators are naturally orthogonal, which can be considered as 
a special case of both unitary and real-similarity transformations. As a consequence, only 
(complex) symmetric Hamiltonians are constructed. Thus, the VT symmetry condition and 
P-pseudo- Hermit icity condition coincide, 

VH = H*V VH = H ] V. (77) 

Although this limits the number of parameters in the Hamiltonian matrices, it has the 
advantage to have both the VT eigenstates and V being an indefinite metric operator. To 
find the proper metric W, one may construct the so-called C operator which commutes with 



both H and VT 17J]. By construction, the eigenvalues of C are set to cancel the negative 
signs in the P-norm. In this case, the metric operator can be chosen as W = VC. 

If we look at the relation between two formulations from a different angle: For a VT- 
symmetric matrix Hamiltonian H, can we find a V operator (V ^ V in general) such that 
H is also "P-pseudo-Hermitian? Or for a P-pseudo-Hermitian matrix Hamiltonian H, can 
we find a parity V such that H is also PT-symmetric? These questions are much harder to 
answer. 

To answer these questions, let us recall a theorem from linear algebra: For any given 
n x n matrix B, there exists a similarity transformation, Ab, which transposes B. That 
is, AbBA^ 1 = B T . A simple proof of this theorem can be found in|A] Equipped with this 
theorem, it is straight forward to verify that a "PT-symmetric H must satisfy 

QH = H ] Q with Q = A* H V, (78) 
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where Ah is the transpose matrix of H: AhHAJ^ = H T . For any given H, there are many 
Ahs which can transpose it. Among these Ahs, if some Ahs make Q Hermitian, Q = Q\ 
then H is Q-pseudo-Hermitian. If one can further find an Ah such that Q is not only 
Hermitian but also a involution, Q 2 = 1, then Q can be considered as the V operator in 
P-pseudo-Hermiticity discussed in this section. We will show that in the case of 2 x 2, such 
a P operator always exists. That is, any 2x2 PT-symmetric matrix Hamiltonian is always 
P-pseudo-Hermitian for some P. 

The reverse problem is similar. For a P-pseudo-Hermitian matrix Hamiltonian H, it 
satisfies 



If one can find an Ah such that Q = Q* and Q 2 = 1, then Q can be used as the parity 
operator in the previous section. Namely, H is PT-symmetric. Again, in 2 x 2, this always 
can be done. 



In this subsection, we show explicit examples of 2 x 2 pseudo-Hermitian matrices. We 
start with Vo to be the third Pauli Matrix, 



QH = H*Q with Q = (A^y 1 P, 



(79) 



D. 2x2 




(80) 



A Po-pseudo-Hermitian Hamiltonian has the form 




(81) 



The eigenvalues of Hq are 



E ± = e ± vV - P 2 - 



(82) 



They are real when r ) 2 > p 2 . The eigenvectors are 




(83) 
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where n± is the normalization constant. Solving the self-adjoint condition in Eq. ([T]), we get 
the form of the metric operator, 

(7 + v pe [S \ 
, (84) 
pe~ lS 7 — v I 

where u and v are arbitrary constants with the constraints w-f > and v 2 < 7 2 — p 2 . The 
eigenvalues of Wo are 

ou± = u (7 ± Vp 2 + • (85) 

They are positive definite. With this metric operator, the eigenvectors of H with different 
eigenvalues are orthogonal 

(E+\W \E-) = = (E-\W \E+) (86) 

and they are normalized to 

A4 = (E±\W \E ± ) = 2\n ± \ 2 u^ff^ ± ^ff^) {yi^ 2 ±v) . (87) 

Given u and v, one can always normalize the eigenvector by choosing n± properly. It is also 
true in the reversed way, for arbitrary non- vanishing n±, one can always tune u and v in the 
metric to normalized the eigenvectors. 

A generic unitary matrix in the coset G\(C 2 ) can be parameterized as 

U(9,(p) = e -%"n/* e -V°*/2. (88) 

It transforms Vo to 

- . / cos9 e^sinA 

V = U{9, p)V U\9, <p) = n r ■ a = , (89) 

V e lv sin# -cos9 / 

where n r = (sin# cosp, sin# sin tp, cos 8) is a unit vector. The P-pseudo-Hermitian Hamil- 
tonian can be obtained by the same unitary transformation, 

H = U(9,<p)H rf(9,<p) = ea + (7 n r + ip sin 6 n 9 + ip cos 6 n v ) • <x, (90) 

where n" = (cos 9 cos p>, cos 9 sin ip, — sin 9) and = (— simp, cos ip, 0) are two unit vectors 
perpendicular to n r . Finally, the metric operator for this H is 

W = U{9,p)WoU\9,p) =u [ 7 a + (vn r + pcosS n e — p sin 5 n^) ■ cr] . (91) 
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Now let us see what goes wrong when we approach a Jordan block. Unitary transforma- 
tions do not change the diagonalizability of a matrix. For simplicity, we work in frame. 
H forms Jordan block if and only if 7 2 = p 2 ^ 0. Define a small positive parameter e by 

P 2 =7 2 (l-e). (92) 

To avoid double-limit, we set v = 0. For a fixed u, when e is small, the larger eigenvalue of 
Wo, u > behave well but the smaller one, u < is vanishing, 

u > ~ 2^7, and u < ~ \u^e. (93) 

A similar ill behavior can be found in the normalization of eigenvectors of H , 

A4 ~ 2|n±| 2 M7 3 e. (94) 

When Hq is a Jordan block, one can only find one eigenvector. Another vector lives in 
the Jordan chain. For simplicity, let us assume 7 = p 7^ 0. In this case, H has the form 

( e + 7 7e i<5 \ , , 

H =[ 11 . (95) 

\ — 7e e — 7 / 
The eigenvalue of Hq is e and the only eigenvector is 

|$o)=n ^ ^j, (96) 

where no is the normalization constant. The second vector in the Jordan chain satisfies 

(H -ea Q )\$ 1 ) = \® Q ). (97) 

The solution is easy to find, 

\®i) = n ( ° +a|$ ), (98) 

with a to be an arbitrary constant. 

All 2x2 PT-symmetric matrix Hamiltonians are also pseudo-Hermitian with respect 
to some V operators and vice versa. For example, the Po7~- s y mme tric H in Eq. (126]) is 
P-pseudo-Hermitian with 

(cos 5 — i sin 5 , 
I • (99) 
i sin 5 — cos 5 
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The Po-pseudo-Hermitian Hq in Eq. ([Mil is also PT-symmetric with 

1 / 7 pcos5 . 

v = / 2 2 2 A / i- 



One more example, the parity operator considered in Ref. 
Eq. ([36]) with ip = 0: 



131 ] is the special case of Pi in 



■ cos 9 sin 9 

Vi = Ri(e,o)v R[ \e,o)= ( |. (101; 

sin — cos 

The PxT-symmetric matrix takes the form of Hi in Eq. ( 1381) with p = 0. This f/i is 
P-pseudo-Hermitian with 

(cos 5 cos # cos 5 sin # — i sin 5 \ 
(102) 
cos S sin # + i sin 5 — cos S cos # / 

In Appendix [Bj we give the explicit form of V for a generic 2x2 PT-symmetric matrices in 
Eq. (138|) and Eq. (|39l) to be P-pseudo-Hermitian. We also give the explicit form of the parity 
operator V for a generic 2x2 P-pseudo-Hermitian matrices in Eq. (|90|) to be PT-symmetric. 



IV. GENERALIZED VT SYMMETRY 

We generalize PT symmetry by defining only the combined PT operation as an anti- 
unitary involutionary operator: 

VT = V*, and (PT) 2 = 1, (103) 

where * represents the complex conjugation just as in Sec. [Til These conditions immediately 
leads to 

VP* = t. (104) 

Comparing to PT symmetry, the generalized PT symmetry lifts the restriction of the parity 
operator to be an involution as well as removes the commuting relation between the parity 
and the time reversal operators. Meanwhile, we maintain the combined operation to be an 
involution. As shown in Sec. [Uj the combined involutory condition and the anti-linear prop- 
erty are what essentially lead to a possible real spectrum and the phenomena of spontaneous 
PT symmetry breaking. A generalized PT-symmetric matrix Hamiltonian satisfies 

[H,VT} = VH* = HV. (105) 
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The transformations on the V matrix that preserves Eq. (11041) are 



V = AP (A^)\ 



(106) 



Under this transformation, a generalized "PoT-symmetric matrix Hamiltonian Hq transforms 
to a generalized PT-symmetric matrix Hamiltonian H if 



which turns out to be a similarity transformation. There is no limitation on the transfor- 
mation matrix A except it must be invertible. For the case of iV x N, A e GL(JV, C). 

When the eigenvalues of H are real, the argument in Eq. ffHJj) still holds. That is, the 
eigenstates of H are also eigenstates of VT with the eigenvalues can be chosen to be one. 

If Hq not only has an entire real spectrum, but also is diagonalizable, one may solve the 
self-adjoint condition in Eq. (Tj[|) to find the proper metric operator Wq. The corresponding 
metric for H can be obtained by 



The similarity transformation in Eq. (j!07p also means that all self-adjoint matrices are 
generalized PT-symmetric. There are two ways to see it. The simpler way is to start with a 
diagonal form of Hq. It is obviously generalized PoT-symmetric with Vq = 1. The similarity 
transformation brings H to be an arbitrary self-adjoint matrices with the same spectrum. 
This H is generalized PT-symmetric with P = A (A -1 )* . Note that although Vq is the unity, 
V is non-trivial in general: V ^ 1. 

Another way is to start with a diagonal metric operator Wd' 



H = AiJnA" 1 



(107) 



w = (A-ywoA -1 . 



(108) 



W D = Diag{wi, w 2 , w 3> • ■ •}• 



(109) 



The self-adjoint condition in Eq. (P) leads to the Hamiltonian with the form 



/ 




H D 



(ai2-i6i2) a 22 Z3 J^-(a 23 + i& 23 ) ••• 

(ai 3 -i&i 3 ) ^S-(a23-i&2 3 ) ^33 



(110) 



V 
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with all djj and bij to be real. We observe three properties for a matrix to be self-adjoint 
respect to a diagonal metric: (i) All the diagonal elements must be real, 



H 



D 



( h d)*\ 



111) 



(ii) The off-diagonal elements have opposite phases comparing to the mirror element on the 
other side of the diagonal, 

h d h d = ( h d h dY ! (H2) 

(iii) And the ratio between two moduli of mirroring off-diagonal elements is determined by 
the eigenvalues of Wd, 



TTl] 

n D 



TT'Jt 

D 



ijj, 



(113) 



Note that Hd is a Hermitian matrix if and only if Wd proportional to the unity matrix. 

Now let us see what a generalized PoT-symmetric matrix with a diagonal Vq looks like. 
Because of Eq. (11041) . a diagonal Vq has matrix elements to be a pure phase on the diagonal: 



V = Diag{e ia \e iQ2 ,e ia V--}. 



114) 



Solving the generalized VT symmetry condition in Eq. HI 05 j) . we get the matrix elements of 
H satisfies 



Multiply e 1 ( a i+ a i)/ 2 on the both sides, we have 



(115) 



(116) 



This means that the phase of the ij element of Hq is (at — atj)/2. Thus Hq has the form 

/ t-,, ri2e i("i- a 2)/2 ri3e i(ai-a 3 )/2 ...\ 



Til 



i(a 2 -ai)/2 roop i(a 2 -a 3 )/2 

^33 



r3ie i(aa-ai)/2 



?"22 r 23 e J 
i(a3-a2)/2 



V 



117) 



with Tij real. The form of the generalized PoT-symmetric H is very similar to the self- 
adjoint Hd in Eq. (II 101) . But H in Eq. (11171) is more general. The subtle difference lies in 
the third property of Hd in Eq. f II 1 3 j) : Hp and H D % must be both zero or both non- vanishing. 
On the contrary, one can take Hq = at the same time Hq ^ for i ^ j. In addition, 
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the eigenvalues of Hp must be real and Hq allows both real and complex conjugate pair 
eigenvalues. We conclude that H is more general than Hd, and Hp is a special case of 
H . Finally, the transformations in Eq. ( I107p gives us that all self-adjoint matrices process 
a generalized VT symmetry. 

It is very easy to count how many real parameters in an arbitrary N x N self-adjoint 
matrix. A generic N x N complex matrix has 2N 2 real parameters. The reality of the 
eigenvalue poses N constraints. So an iVx iV self-adjoint matrix has 2N 2 —N real parameters. 

Although all self-adjoint matrices form a subset of generalized PT-symmetric matrices, 
the latter have the same number of parameters. It is simply due to the anti-linear involutory 
property of VT, a generalized PT-symmetric matrix have real or complex pair eigenvalues. 
For an Af x iV matrix, it poses N constraints. 

Reducing the generalized VT symmetry to the VT symmetry discussed in Sec. [Til is 
trivial, all we need is to start with a real Vq = Vo and transform only by real matrices. 
In this way, the resulting parity operator V is always real and H is PT-symmetric. On 
the contrary, reducing the generalized VT symmetry to pseudo-Hermiticity is much more 
complicated. For a generalized "PT-symmetric matrix H, VH* = HV, if there exist a matrix 
Ah such that 

A H HA H l = H T and A H A* H = 1, (118) 

then H satisfies 

QH ] = HQ with Q = VA H . (119) 

Among these A^s, if some make Q Hermitian, Q = , then H is Q-pseudo-Hermitian. 
Finally, if one can find an Ah such that Q 2 = 1 also holds, then this Q matrix can be the 
V operator discussed in Sec. that is, H is "P-pseudo-Hermitian. 

For completeness, we give an explicit parametrization of the 2 x 2 V operator discussed 
in this section: 

(cos 9 + i sin 9 sin 5 i (sin 9 cos 5 — tanh ip) 
K ' | . (120) 

i (sin 9 cos 5 + tanh <p) cos 9 — i sin 9 sin 5 

It is straight forward to verify that VP* = cr . 
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TABLE I: The numbers of real parameters in real-symmetric, Hermitian, "PT-symmetric, P-pseudo- 
Hermitian, self-adjoint, and generalized "PT-symmetric matrices. 



Dimension 


2 3 m + n 


A^> 1 


Real-symmetric 


3 6 n(m + n)(m + n + 1) 


iiV 2 


Hermitian 


4 9 (m + n) 2 


N 2 


PT-symmetric or 


6 13 (m + re) 2 + 2mn 


N 2 ~ |A^ 2 


P-pseudo-Hermitian 




Self-adjoint or 


6 15 2(m + n) 2 — (m + n) 


2iV 2 


Generalized PT-symmetric 







V. CONCLUSIONS 

In conclusion, we list the parameter counting of all relevant matrices in Table [H From the 
table, in the case of 2x2, PT-symmetric matrices, P-pseudo-Hermitian matrices, generalized 
■PT-symmetric matrices, and self-adjoint matrices all have the same number of parameters. 
In fact, any 2x2 self-adjoint matrix is PT-symmetric, "P-pseudo-Hermitian, and generalized 
"PT-symmetric. This is not true in another way around because VT symmetry may be 
broken. Even it does not, a PT-symmetric matrix or a P'-pseudo- Hermitian or a generalized 
"PT-symmetric matrix could form a Jordan block and no longer diagonalizable. In the 
dimension higher than 2, a self-adjoint matrices is always generalized "PT-symmetric but 
not PT-symmetric or "P-pseudo-Hermitian in general. 

Appendix A: Proof of the existence of the transpose matrix 

Suppose that an re x re matrix B is similar to a Jordan block J n , 

FBF- 1 = J n . (Al) 
Consider another re x re matrix As defined as 

A B = F T S n F, (A2) 



24 



where S n is an SIP defined in Eq. ( )67|) . then Ab is the transpose matrix of B, 

A B BA B l = B T . (A3) 

If matrix B is similar to a direct sum of Jordan blocks, 

FBF" 1 = J mi © J ni2 © • • • , (A4) 

then the transpose matrix can be constructed by replacing S n in Eq. flA2j) by a similar sum, 
S mi © S m2 © • • •. Since any square matrix is either similar to a Jordan block or a direct sum 
of Jordan blocks, there always exists a transpose matrix as constructed above. 



Appendix B 

In this appendix, we give the explicit forms in 2 x 2 about the equivalence of VT symmetry 
and V pseudo-Hermiticity. 

The generic "P-pseudo-Hermitian H in Eq. fl§D|) is PT-symmetric respect to 



y V21 —7 cos 9 
with 

"P 12 = 7 sin # cos <p + p cos 5 cos <p + p sin 5 cos # sin <p 
V21 = 7 sin 6* cos <p — p cos 5 cos </? — p sin 5 cos # sin tp 

Ai = 7 2 (cos 5 cos # sin <p> — sin 5 cos ip>) 2 + (7 2 — p 2 ) (sin S cos # sin ip + cos 5 cos ip) 2 . 

The generic "PxT-symmetric Hi in Eq. (|38|) is "P-pseudo-Hermitian respect to the Hermi- 
tian involution operator 

1 / 7 cos 5 cos 9 V12 
V = ±-=- ' _ (B2) 
vA 2 \ p 21 — 7 cos 5 cos # 



with 



V12 = 7 cos 5 sin 9 cosh <p — i7 sin 5 cosh <p — ip cos # sinh <p 
■p 2 i = 7 cos 5 sin 9 cosh <p + i7 sin 5 cosh <p + ip cos 9 sinh tp 
A 2 = 7 2 cos 2 5 (sin 2 # cosh 2 <p + cos 2 9) + (7 sin 5 cosh <p + p cos 9 sinh <p) 2 . 
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The generic T^T-symmetric H2 in Eq. (|39|) is "P-pseudo-Hermitian with respect to the 
Hermitian involution operator 

II 7 cos 5 cosh 9 V[n 
V = ~ (B3) 

v A3 \ p' — 1 cos 8 cosh 6* 



with 



7^2 — ~~ 7 cos 5 sinh 6* sinh <p — i7 sin 5 cosh # cosh <p — ipsinhy? 
V21 = — 7 cos 5 sinh 6* sinh <p + i7 sin 5 cosh # cosh <p + ipsinh(p 
A3 = 7 2 cos 2 <5(l + sinh 2 9 cosh 2 <p) + (7 sin 5 cosh 6* cosh ip + p sinh ip)" 
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